Colored Non-Crossing Euclidean Steiner Forest

Given a set of $k$-colored points in the plane, we consider the problem of finding $k$ trees such that each tree connects all points of one color class, no two trees cross, and the total edge length of the trees is minimized. For $k=1$, this is the well-known Euclidean Steiner tree problem. For general $k$, a $k\rho$-approximation algorithm is known, where $\rho \le 1.21$ is the Steiner ratio. We present a PTAS for $k=2$, a $(5/3+\varepsilon)$-approximation algorithm for $k=3$, and two approximation algorithms for general~$k$, with ratios $O(\sqrt n \log k)$ and $k+\varepsilon$

Multi-messenger observations of a binary neutron star merger

On 2017 August 17 a binary neutron star coalescence candidate (later designated GW170817) with merger time 12:41:04 UTC was observed through gravitational waves by the Advanced LIGO and Advanced Virgo detectors. The Fermi Gamma-ray Burst Monitor independently detected a gamma-ray burst (GRB 170817A) with a time delay of ~1.7 s with respect to the merger time. From the gravitational-wave signal, the source was initially localized to a sky region of 31 deg2 at a luminosity distance of 40+8-8 Mpc and with component masses consistent with neutron stars. The component masses were later measured to be in the range 0.86 to 2.26 Mo. An extensive observing campaign was launched across the electromagnetic spectrum leading to the discovery of a bright optical transient (SSS17a, now with the IAU identification of AT 2017gfo) in NGC 4993 (at ~40 Mpc) less than 11 hours after the merger by the One- Meter, Two Hemisphere (1M2H) team using the 1 m Swope Telescope. The optical transient was independently detected by multiple teams within an hour. Subsequent observations targeted the object and its environment. Early ultraviolet observations revealed a blue transient that faded within 48 hours. Optical and infrared observations showed a redward evolution over ~10 days. Following early non-detections, X-ray and radio emission were discovered at the transient’s position ~9 and ~16 days, respectively, after the merger. Both the X-ray and radio emission likely arise from a physical process that is distinct from the one that generates the UV/optical/near-infrared emission. No ultra-high-energy gamma-rays and no neutrino candidates consistent with the source were found in follow-up searches. These observations support the hypothesis that GW170817 was produced by the merger of two neutron stars in NGC4993 followed by a short gamma-ray burst (GRB 170817A) and a kilonova/macronova powered by the radioactive decay of r-process nuclei synthesized in the ejecta

Brazil nut's effect beyond spherical grains: Elongation matters!

This is the first experimental study of segregation of mixtures of spherocylindrical rods of same volume and different elongations (down to the limit of spheres, where elongation is zero). When vibrated in a cylinder, we observe that short spherocylinders tend to migrate above isovolume spheres, while long spherocylinders tend to migrate to the bottom. More generally, short spherocylinders tend to migrate to the top of long spherocylinders. Moreover, the particles migrating to the bottom tend to arrange in an ordered state. We couple experiments and numerical simulations. Our discrete element model (DEM) faithfully captures the physical behavior, contrarily to previous studies made with Monte Carlo methods. We extend the study to binary mixtures of spheres and spherocylinders that have different volumes. We link these observations to the tendency of rods to orient vertically and interpret both phenomena through an extension of the void-filling mechanism

Brazil nut's effect beyond spherical grains: Elongation matters!

This is the first experimental study of segregation of mixtures of spherocylindrical rods of same volume and different elongations (down to the limit of spheres, where elongation is zero). When vibrated in a cylinder, we observe that short spherocylinders tend to migrate above isovolume spheres, while long spherocylinders tend to migrate to the bottom. More generally, short spherocylinders tend to migrate to the top of long spherocylinders. Moreover, the particles migrating to the bottom tend to arrange in an ordered state. We couple experiments and numerical simulations. Our discrete element model (DEM) faithfully captures the physical behavior, contrarily to previous studies made with Monte Carlo methods. We extend the study to binary mixtures of spheres and spherocylinders that have different volumes. We link these observations to the tendency of rods to orient vertically and interpret both phenomena through an extension of the void-filling mechanism

Vertical ordering of rods under vertical vibration.

Granular media composed of elongated particles rearrange and order vertically upon vertical vibration. We perform pseudo-two-dimensional discrete element model simulations and show that this phenomenon also takes place with no help from vertical walls. We quantitatively analyze the sizes of voids forming during vibrations and consider a void-filling mechanism to explain the observed vertical ordering. Void filling can explain why short rods are less prone to align vertically than long ones. We cannot, however, explain, invoking just void filling, the existence of an optimum acceleration to promote vertical ordering and its dependence on particle length. We finally introduce an interpretation of the phenomenon, by considering the energetic barriers that particles have to overcome to exit a horizontal or a vertical lattice. By comparing these energetic thresholds with the peak mean particle fluctuant kinetic energy, we identify three different regimes. In the intermediate regime a vertical lattice is stable, while a horizontal is not. This interpretation succeeds in reconciling both dependencies on vibration acceleration and on particle length

About the influence of friction and polydispersity on the jamming behavior of bead assemblies.

We study the jamming of bead assemblies placed in a cylindrical container whose bottom is pierced with a circular hole. Their jamming behavior is quantified here by the median jamming diameter, that is the diameter of the hole for which the jamming probability is 0.5. Median jamming diameters of monodisperse assemblies are obtained numerically using the Distinct Element Method and experimentally with steel beads. We obtain good agreement between numerical and experimental results. The influence of friction is then investigated. In particular, the formation of concentric bead rings is observed for low frictions. We identify this phenomenon as a boundary effect and study its influence on jamming. Relying on measures obtained from simulation runs, the median jamming diameter of bidisperse bead assemblies is finally found to depend only on the volume-average diameter of their constituting beads. We formulate this as a tentative law and validate it using bidisperse assemblies of steel beads

On non-rank facets of the stable set polytope of claw-free graphs and circulant graphs

We deal with non-rank facets of the stable set polytope of claw-free graphs. We extend results of Giles and Trotter [7] by (i) showing that for any nonnegative integer a there exists a circulant graph whose stable set polytope has a facet-inducing inequality with (a,a+1)-valued coefficients (rank facets have only coefficients 0, 1), and (ii) providing new facets of the stable set polytope with up to five different non-zero coefficients for claw-free graphs. We prove that coefficients have to be consecutive in any facet with exactly two different non-zero coefficients (assuming they are relatively prime). Last but not least, we present a complete description of the stable set polytope for graphs with stability number 2, already observed by Cook [3] and Shepherd [18]

Coupled tomography and distinct-element-method approach to exploring the granular media microstructure in a jamming hourglass.

We describe an approach for exploring microscopic properties of granular media that couples x-ray microtomography and distinct-element-method (DEM) simulations through image analysis. We illustrate it via the study of the intriguing phenomenon of instant arching in an hourglass (in our case a cylinder filled with a polydisperse mixture of glass beads that has a small circular shutter in the bottom). X-ray tomography provides three-dimensional snapshots of the microscopic conditions of the system both prior to opening the shutter, and thereafter, once jamming is completed. The process time in between is bridged using DEM simulation, which settles to positions in remarkably good agreement with the x-ray images. Specifically designed image analysis procedures accurately extract the geometrical information, i.e., the positions and sizes of the beads, from the raw x-ray tomographs, and compress the data representation from initially 5 gigabytes to a few tens of kilobytes per tomograph. The scope of the approach is explored through a sensitivity analysis to input data perturbations in both bead sizes and positions. We establish that accuracy of size--much more than position--estimates is critical, thus explaining the difficulty in considering a mixture of beads of different sizes. We further point to limits in the replication ability of granular flows away from equilibrium; i.e., the difficulty of numerically reproducing chaotic motion

Coupled tomography and distinct-element-method approach to exploring the granular media microstructure in a jamming hourglass.

We describe an approach for exploring microscopic properties of granular media that couples x-ray microtomography and distinct-element-method (DEM) simulations through image analysis. We illustrate it via the study of the intriguing phenomenon of instant arching in an hourglass (in our case a cylinder filled with a polydisperse mixture of glass beads that has a small circular shutter in the bottom). X-ray tomography provides three-dimensional snapshots of the microscopic conditions of the system both prior to opening the shutter, and thereafter, once jamming is completed. The process time in between is bridged using DEM simulation, which settles to positions in remarkably good agreement with the x-ray images. Specifically designed image analysis procedures accurately extract the geometrical information, i.e., the positions and sizes of the beads, from the raw x-ray tomographs, and compress the data representation from initially 5 gigabytes to a few tens of kilobytes per tomograph. The scope of the approach is explored through a sensitivity analysis to input data perturbations in both bead sizes and positions. We establish that accuracy of size--much more than position--estimates is critical, thus explaining the difficulty in considering a mixture of beads of different sizes. We further point to limits in the replication ability of granular flows away from equilibrium; i.e., the difficulty of numerically reproducing chaotic motion